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1 Introduction 



The problem of constructing 3-manifolds of prescribed non-negative scalar curvature is im- 
portant for the initial value formulation of the Einstein equations in general relativity in the 
asymptotically flat case. For some standard references see (2j[Hl|6]. For a more recent survey 
article involving the constraint equations see [3]. In the maximal "gauge", valid initial data 
consist of an asymptotically flat Riemannian manifold (M, g) together with a symmetric, 
trace free tensorfield kij, vectorfield J, and non-negative function p. Here k^ is interpreted 
as the second fundamental form of (M, g) as embedded in spacetime, J is the mass current 
density, and p is the local mass density; the non-negativity of the latter follows from the 
dominant energy condition. In addition, the data (M, g, k, J, p) must satisfy the Einstein 
constraint equations, which in the maximal "gauge" appear as 



where R(g) denotes the scalar curvature of the metric g and V refers to the Levi-Civita 
connection with respect to g. Thus, the construction of such initial data involves the con- 
struction of a manifold (M, g) of non-negative scalar curvature. In the time symmetric case 
k = the scalar curvature is proportional to the mass density R{g) = lQirp, and the problem 
reduces to a prescribed non-negative scalar curvature problem. 

The prevalent method for solving the constraints has been the conformal method. This 
requires at the outset a metric g which is conformally equivalent to a metric of non-negative 
scalar curvature. Such a metric g is only guaranteed to exist if the L 3//2 norm of the negative 
part of the scalar curvature R(g) is small [TJ. That is, the scalar curvature of g must be almost 
non-negative already. Thus, what is needed is a non-conformal method for constructing 
general metrics of prescribed non-negative scalar curvature. Solving the parabolic scalar 
curvature equation (II. 4p below provides just such a method. It is our main goal, in the 
present paper, to construct many solutions of the parabolic curvature equation by means of 
equivariant bifurcation theory and symmetry breaking. For previous ad-hoc constructions 
see [32] . Figure 11.11 sketches a 2-dimensional caricature of the 3-dimensional space initial 
data for a black hole M foliated over a radial variable r by 2-dimensional spheres S = S 2 , 
which are only horizontal circle in Figure 11.11 

Given a foliation of the manifold M, with foliating function r, the parabolic scalar curvature 
equation relates the scalar curvature to the transversal metric component u = \Vgr\~ 1 . We 
will only be concerned with regions upon which r is non-critical so that we may assume 
without loss of generality that the manifold takes the form M = [ro,ri] x S, where £ is a 
regular 2-manifold - in the present work we take £ = S 2 , but the equation is true for any E. 
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Figure 1.1: 2d- caricature of a black hole initial data for the Einstein equations. The vertical 
variable is the radius r. Each horizontal section {r = constant} in this figure is a circle, 
standing in fact for a 2-sphere E = S 2 . The region outside the apparent horizon (and above 
it in this picture ) was constructed in f3^ \3J$ . The present paper constructs the region called 
M between an initial radius ro > and a blow up radius r\ = 1 whose section {r±} x S 2 is a 
critical point of the area functional. The matching and the center regions (shaded) are still 
to be constructed. 



After a foliation preserving diffeomorphism, any metric g on M can be expressed as 
(1.3) g = u 2 dr 2 + r 2 u, 

where u is a family of metrics on E extended to M so that u(d r , •) = 0. The parabolic scalar 
curvature equation can then be expressed as 



(1.4) Hrd r u = u 2 A w u + Bu - - -r 2 Rj u\ 
where 

(1.5) H = 2 + rd r u ab u a \ 

A w is the Laplace-Beltrami operator of u, k is the Gauss curvature of u, and B is also 
determined entirely by u; specifically B = rd r H — H + \{H 2 + \u + rd^l^). Equation fll.4p 



was first derived in the quasi-spherical case by Bartnik [3J. It was fully generalized to the 
3-dimensional case by Smith and Weinstein [35]. It was later extended to n-dimensions by 
Shi and Tarn [29] . 

Concerning the type of equation (jl.4p . note that the mean curvature of the foliation is given 
by H = H/ru. Thus, assuming that the family u is such that the mean curvature of the 
foliation is positive, equation fll.4p is in fact parabolic, with the radial variable r playing the 
role of the time variable. Accordingly, we shall often speak of r and increasing functions of 
r as time variables in analogy with the heat equation, porous medium equation, etc. For 
convenience, we also call the boundary value problem with prescribed data at {ro} x E an 
initial value problem. 

Assuming compactness of E, one uses this equation to construct a metric g of prescribed 
scalar curvature R in the following way: choose a smooth family of metrics such that H > 0, 
as well as "initial" data uq G C°°(E), < uq. Solve the parabolic scalar curvature equation 
f ll.4j) with initial data u(tq) = Uq. By standard parabolic regularity theory this is always 
possible for a small enough interval [ro, tq + e). 

Of course, even in the case of long time existence, the resulting manifold will not be complete 
since we have a boundary at r = r . Assuming E = S 2 , it is possible to produce solutions 
starting at r = [3J, but we shall not worry about this in the present work. For small 
r G [0, ro] it is simpler for the metric to have been constructed by other means. For instance, 
one can first construct the metric in Gaussian normal coordinates for r < ro, and then use 
the parabolic scalar curvature equation to extend to r > r . For more interesting topologies, 
one will need to start with an inner minimal surface, which will then be joined to a similarly 
constructed manifold on the other side. For more on this see [33] , [31] . 
Concerning global existence, it suffices that the operator 

(1.6) T = A u -L- r ^\ 

be non-positive [33J . Since we are primarily interested in asymptotically flat manifolds whose 
ends have the topological structure of [r , +oo) x S 2 , we choose, in fact, E = S 2 . Then by 
the Gauss-Bonnet theorem J s2 k = An, and we are allowed to choose to such that k > 0. 
Positivity of the operator — T can be ensured by choosing R < 2r 2 K. 

Nonetheless, there are many situations in which this condition is not satisfied that do, in 
fact, lead to blow-up at finite r. If the mass density p = R/16ir is very large so that, for 
instance k— < — c < 0, one can easily see from the maximum principle that blow-up will 
occur at finite time (radius) r%. This is, of course, only blow-up of a metric component, but 
it is geometrically significant with our choice of foliation: the mean curvature of the foliation 
approaches at any blow-up point, and parabolicity is no longer satisfied at the blow-up 
"time" r = t\. This will be true for any bounded "time" variable that we might use. 
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It may, however, happen that the constructed metric itself retains sufficient regularity at the 
blow-up "time" r = r% so as to allow an extension beyond this radius by other means. The 
simplest example occurs when u is just the standard metric u of S 2 and r 2 R = 4. Then the 
parabolic scalar curvature equation becomes 

(1.7) 2rd r u = u 2 Au + u 3 + u, 

where we have written A = On [r , 1) x S 2 , for instance, this admits the trivial blow-up 
solution 

(1.8) u = (l/r-l)- 1/2 , 
with corresponding metric 

(1.9) g= (1/r- l)~ 1 dr 2 + r 2 Cj 

Obviously the metric blows up at r = 1. In terms of the geodesic distance function 

(1.10) s = f (1/r - l)~ 1/2 dr~, 

Jr 

however the metric takes the form g = ds 2 + r 2 (s)u. It is easily checked that r is C°° smooth 
in s, and thus g is C°° smooth up to and including the boundary at r = 1. 
A natural next line of investigation, which we continue to follow in this work, is to ask if 
there exist anisotropic solutions u of (11.41) in the case u = u with the same property. Here 
anisotropic means, in other words, that u = u(r,p) is not only a function of r alone but 
also depends on the angular variables p G X = S 2 . Such nonhomogeneous solutions are 
not invariant under the action of the full symmetry group 0(3) any more: they exhibit 
anisotropy. More specifically, we would like to address the existence of self-similar blow-up 
solutions in the form 

(1.11) u= (1/r- iy 1/2 v, 

for a function v = u(r,p), which is bounded above and below by positive constants, and is 
in addition C°° on the interior (r , 1) x S 2 and C° at the boundary {ro, 1} x S 2 . Using the 
function s defined in the previous paragraph (which is now only proportional to the geodesic 
distance), the corresponding metric g can be expressed as 

(1.12) g = u 2 ds 2 + r 2 (s)u, 

and the aforementioned regularity of v is exactly reflected in the metric g. In this work we 
shall not always obtain regularity at the boundary r = 1 beyond C°, but it should in any 
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case be noted that in order to do so it is more appropriate to use a function such as s for 
the foliating function since the area radius variable r is degenerate at r = 1. In fact, in the 
present case we can say more: the mean curvature H of the outer boundary surface r = 1 
vanishes identically; i.e. it is a minimal surface. Even in the case that further differentiability 
at r = 1 is lacking, the stability functional - alias the second variation of the area functional- 
is defined and shows this surface to have at least one unstable direction - more on this in 
Section 7. 

The approach described in the preceding paragraph was initiated in [32]. In that work it is 
shown that if the function v is bounded away from 0, then it is also uniformly bounded above, 
and thus the final constructed metric is C°°([ro, 1) x S 2 ) D C°([ro,l] x S 2 ). Furthermore, 
under suitable conditions on r 2 R, non-linear stability of the trivial blow-up was obtained. In 
particular this showed the existence of at least some non-trivial solutions of fll.4p whose blow 
up behavior is only asymptotically given by the trivial solution f II . 8 [) . However, this could 
certainly not be described as a plethora of solutions. Part of the aim of the present work is 
to find more blow-up solutions in this category by a more systematic bifurcation analysis. 
Define A = (r 2 R — 2). In this work we take A to be an adjustable constant, i.e. a bifurcation 
parameter. Substituting (II. lip u = (1/r — \)~ x l 2 v in equation ( II. 4p we obtain the equation 
for v: 

(1.13) 2(1 - r)d r v = u 2 Au + \? - v. 

We may eliminate the breakdown of parabolicity at r = 1 by defining a new variable t such 
that r = 1 — exp(— 2t). The surfaces r — > 1 are described by t — > oo in the self-similarly 
rescaled equation. 

(1.14) d t v = v 2 Av + \? - v. 

A few observations are in order: C° Regularity of the metric g at r = 1 corresponds to 
the existence of a globally bounded solution u, < \i < v < M < oo, converging to an 
equilibrium G C°°(S 2 ) of equation (j!.14p as t — > oo. The equilibrium is a solution of 
the equation 

A 1 

(1.15) Az/* + -v* = 0. 

<-< 

It will be convenient to rescale, v t— > vj a/A/2, t t— > t\/2, so that the equilibrium equation 
( I1.15P becomes 

(1.16) Aza 
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and the dynamical equation becomes 



(1.17) d t u = u 2 Uv+~(v-^ 

In order to apply analytic semigroup theory and equivariant bifurcation theory, it is more 
convenient to translate the trivial equilibrium solution z/* = 1 to the origin by defining 
v — v — 1: 




(1.18) 8 t v = (v + l) 2 (Av + Xf(v)) . 
Here 

(1.19) f(v)=v-ff l /(l + v). 

We often refer to equation (jl.l8p as the rescaled equation in the following. The equilibria 
of the rescaled equation verify 

(1.20) Au, + A/(u,) = 0. 

Note that such an equilibrium solution corresponds to self-similar blow-up 

(1.21) u(r,p) = (1/r- iy 1/2 u,(p), peS 2 . 

Rephrased in the new notation, the work [32] only treats the case A = 1 explicitly: trivial 
self-similar blow-up corresponding to v* = 0, and a local strong stable manifold of solutions 
which converge to this. In the present work we vary A and use 0(3) equivariant bifurcation 
theory to produce branches of solutions. Indeed, as A crosses the eigenvalues \i = £(£ + 1) of 
— A52, the linearized left hand side of the equation fll.20p becomes linearly degenerate, and 
0(3) equivariant bifurcation theory produces branches of solutions along the appropriate 
isotropy subspaces. At each branch point we obtain one or several local symmetry breaking 
families of solutions v* of the equilibrium equation. Each such solution v* yields a self-similar 
metric metric 

(1.22) g = l!±0L ir i +r ^ 

which can be seen to be C°° smooth up to and including the boundary r = 1 in the radi- 
ally geodesic coordinates s. These anisotropic metrics do not posses full 0(3) symmetry, 
but rather the symmetries of the isotropy subspaces of the eigenspace of A52, which were 
encountered at the branch point from which it originated. 
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Based on [23], we adapt a strong stable manifold theorem to apply at each of these new 
equilibria. This allows us to construct immortal and eternal solutions v to the translated 
rescaled equation fll . 18j) . Immortal, i.e. defined and bounded uniformly for t — > +00, alias 
r — > 1, are solutions v in the strong stable manifolds of equilibria u*. These converge to 
and correspond to asymptotically self-similar solutions u, as indicated in (11.21) II . 22[) above. 
Eternal solutions v exist for all real t with uniform bounds, and are heteroclinic between 
different equilibria: v — > v± for t — > ±00. The metric interpretation u of ancient solutions v 
which exist and are globally bounded for t — > — 00, alias r — > — 00, and which constitute the 
unstable manifolds of equilibria v of (11.181) . will be discussed in Section 7. 
Equilibria v* and solutions t> asymptotic to them in fact exhaust the dynamic possibilities 
of uniformly bounded eternal solutions to fll. IS)) . Indeed the standard energy functional 

(1.23) E(v):= I ( l\Vv\ 2 - XF(v) J dp, 
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where F is a primitive function of the nonlinearity /, i.e. d v F = f. The functional F is a 
Lyapunov function for the semilinear variant 

(1.24) d t v = Av + \f(v) 

of (11.181) . Indeed (11.241) can be interpreted as the L 2 gradient semi-flow of the Lyapunov 
functional E(y) because 

(1.25) jE{v) = - J^v 2 dp. 

In particular, uniformly bounded solutions tend to equilibrium in any time direction. Sim- 
ilarly, the rescaled curvature equation (II. 18p can be interpreted as the gradient semi-flow of 
E{y) with respect to a slightly adapted L 2 -metric which depends explicitly on v. Specifically 

(1.26) ±E{v) = -J{l + v)- 1 v*dp. 

for solutions v of fll . 18)) with uniformly positive weight (1 + v In particular uniformly 
bounded solutions still tend to equilibrium, in any time direction. See [32] where accumula- 
tion on several equilibria is excluded using |30) . 

Our main results, Proposition 16.2) Theorem 16.11 and Table 17.41 below, can be summarized 
as follows. Equation 11.4) where u is the standard metric on S 2 independently of r and 
A = r 2 R — 2. We obtain blow-up solutions of the form 

_ 1 

(1.27) u(r lP )= ( ^(1/r-l) J * (v(~ log(l - r),p) + lV peS 2 , 



2 w ' V 2 



7 



where v is an anisotropic function on the 2-sphere whose remaining symmetry is described 
by the isotropy groups of Table 16.31 Figure 16.11 shows how those functions bifurcate from 
the fully isotropic constant function on S 2 as the parameter A varies. Table 17.41 describes 
in each case the heteroclinic connections between bifurcating equilibria. The corresponding 
anisotropic metrics on M are regular up to the blow up radius and exhibit there a minimal 
surface. 

In sections 2, 3, we summarize some necessary background from equivariant bifurcation 
theory. In Section 4, we check that one of the main technical ingredients - the strong stable 
manifold - is provided in our quasilinear well as the semi-group framework. In 

Section 5, we specify the isotropy subgroups which play a role in the symmetry breaking 
bifurcations of equation (jl.l8p . The rest of the paper puts all these elements together to 
obtain the aforementioned result; see Sections 6 and 7. 

This work was supported by the Deutsche Forschungsgemeinschaft, SFB 647 "Space-Time- 
Matter" . 

2 Symmetry and equivariance 

We briefly recall some symmetry terminology for equivariant dynamics. In the subsequent 
Section [3], we formulate the equivariant branching lemma which is an elementary but very 
useful variant of a classical bifurcation theorem of Crandall and Rabinowitz [15] . ceFor 
the convenience of the reader the precise role of equivariance will be emphasized. More 
generally see [TSl ESI H2] for a background on symmetry and Lyapunov-Schmidt reduction, 
and [HI E3 ESI E2] for center manifolds in semilinear settings. 

Let X, Y be Banach spaces with bounded linear group actions g x , g Y of the same group T. 
In other words g x , g Y are group homomorphisms from T to the invertible bounded linear 
operators on X, Y, respectively: 

(2.1) £? X (7i-7 2 ) = Q X {li)- Q X { l2 ) 

for all 71, 72 G T, and similarly for Y. We call g x , g Y representations of V on X, Y . Frequently 
we will use the abbreviated notation 7« := g x ('~f)u, for 7 G T and u G X. We call a 
representation g x strongly continuous if (7, v) 1— > £>(7)t> is continuous. For finite-dimensional 
X this is equivalent to continuity of homomorphism g x : T — > GL(X) from T to the general 
linear group on X. 
A map 

(2.2) F:X^Y 
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is called T-equivariant under g x , g Y , if F(g x ( , ~f)v) = g Y ( r y)F(v) holds for all 7 G T, t> G X. 
To simplify notation we simply write this requirement as 



(2.3) Ffrv) = jF(v) 

without too much ambiguity. In the important special case X CY where X is a subspace of 
Y, albeit with a possibly stronger norm, it is particularly convenient to forget the distinction 
between g x and g Y , when g x simply restricts g Y to X. 

The intuitive notion of "symmetry" can be made precise in two slightly different ways. Given 
v G X and a representation g x of T we call 

(2.4) r„ := {7 G r I 1V = v}, 

i.e. the group of elements 7 which fix v, the isotropy of v or the stabilizer of v. For example 
T = T : trivially v — possesses full isotropy. Note that 

(2.5) r 7 „ = 7 r\ l7 - 1 ; 

i.e. elements of the same group orbit possess conjugate isotropy. 
Conversely, given any subgroup K of T, we can consider the closed subspace 

(2.6) Fixx(K) := X K := {v G X | 7 w = v for all 7 G X} 

of K-fixed vectors, which we call the fix space of K. For example K < T v is a subgroup of 
the isotropy r„, for any v G Fix(i^). 

Let F : X — > Y be T-equivariant, as above, and fix any subgroup K < Y. Then F restricts 
to 

(2.7) F : X K -> 

i.e. F maps X-fixed vectors 11 6 I to F-fixed vectors F(-y) G K. Indeed, let 7 G F, u G 
X . Then ^fF{v) = F^jv) = F(v) and hence F(v) G y . In this general way the fix 
spaces Fix(X) capture the idea of an "ansatz" to respect some symmetry property K, quite 
efficiently and generally. 

The isotropy subgroups of a representation g x of T on X form a lattice as follows. Vertices are 
the conjugacy classes of isotropy subgroups. Containment of conjugacy classes of subgroups 
defines the lattice structure. 

A representation g v of T on a real or complex Banach space V is called irreducible if g v 
does not restrict to a representation on any nontrivial closed subspace {0} 7^ V 7^ V. In 
other words g v {V) is never contained in V, for any such V. By Schur's Lemma, the only 
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T-equivariant complex linear maps on a finite-dimensional complex irreducible representation 
are the complex multiples of identity. Therefore irreducible representations play a central 
role in spectral analysis. 

For example consider the Sobolev space X = H 2 (S 2 ), Y = L 2 (S 2 ) and the self adjoint 
Laplace-Beltrami operator —A on Y with domain X>(A) = X. Then —A : X — > Y is 
equivariant under the action of the compact group of orthogonal matrices 7 G T := 0(3) on 
v G Y given by 

(2.8) (jv)(x) = (g ¥ (l)v)(x) := v^x). 

The representation restricts to X = H 2 C L 2 . The spectrum spec(— A) is given by real 
eigenvalues < Ao < Ai < . . . with 

(2.9) Xt = £(£ + 1). 

The associated eigenspaces Vi of dimension 2£ + 1 enumerate the irreducible continuous 

representations of the special orthogonal group S0(3), for £ — 0, 1, 2, They are spanned, 

for each £, by the real and imaginary parts of the spherical harmonics 

(2.10) Y tm (0, if) = P™(costf) exp(imtf) 

with m = 0, . . . , £ . The real Legendre functions Pp are essentially m-th derivatives of 
Legendre polynomials Pe(x) with a prefactor (1 — x 2 ) m ^ 2 . The spherical harmonics Yi m are 
normalized to form a complete orthonormal basis of Y — L 2 (S 2 ), in polar or zenith angle 
coordinates < d < it and azimuth < <p < 2ir. Because the Legendre polynomials Pe(x) 
are even/odd for even/odd £, the reflection 7 = — id ^ 5*0(3) acts as multiplication by (—l) e 
on Vg. 

As a first example we claim 

(2.11) Fix^(50(2))= span Y m . 

Here elements of S'0(2), and the rotations around the polar axis $ = of S 2 , shift the 
azimuth (p but keep the polar angle $ fixed. By Fourier expansion with respect to <p this 
proves (12.111) . We call functions v($,(p) = v(-&) in Fix(S'0(2)) axisymmetric. 
Note how the precise isotropy T v < SO(3) of axisymmetric functions v depends on £. Let 
Z c 2 := {±id} denote the center subgroup of T = 0(3) = SO(3) U (-50(3)). Then 

0(2) ©Z!^ for £ > even 
0(2)~ for £>0 odd 
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for non zero v = Y m . The group 0(2)~ refers to an isotropy subgroup such that r„ is 
isomorphic to 0(2) < SO (3) but T v n SO (3) = SO (2) is a subgroup of index 2. Specifically 
0(2)~ are the planar rotations SO (2) in the equatorial plane & = ir/2 around the polar axis, 
and the in-plane reflections <p t— > tpo — <fi via reflections at perpendicular planes. The group 
0(2) < S0(3), in contrast, achieves reflections in the equatorial plane by 180° rotations 
around in-plane axes with •§ = 7r/2. By 0(2) ©Z§ we denote the direct product of 0(2) with 
±id, i.e. 0(2) © Z c 2 = 0(2) U (-0(2)). 

In Figure [2~T1 we sketch the lattices of (conjugacy classes of) isotropy subgroups I\, for the 
representation (12. 8p of T = 0(3) on Ve, £ = 0, . . . , 4. In parentheses we indicate diniMFix(rt,) 
of r„. A complete classification had been initiated by [2T](1984); see also the detailed and 
corrected accounts in the monumental work of Lauterbach [23] as contained in Chossat 
Lauterbach Melbourne [EI] (1991); see also [2"4] . 

A maximal isotropy subgroup T v is a strict subgroup T v ^ T = 0(3) which is maximal in 
the lattice with that property. In the cases I = 1, . . . , 4 of figure 2.1 these coincide with the 
isotropy subgroups T v of real one- dimensional fix spaces, 

(2.13) dim Fix (r„) = 1. 

Obviously condition fl 2 . 1 3 [) implies that the isotropy subgroup is maximal: one-dimensional 
linear representations g of groups admit only actions ^(7) = ±1. Conversely, however, it is 
not true that fix spaces of maximal isotropy subgroups are one-dimensional. The counterex- 
ample with lowest i is i = 12, the octaedral group r„ = O © Zj, and dimFix(r i) ) = 2. 
For further details on the nomenclature of subgroups of 0(3) we refer to [EH H21 EI] • We only 
describe the maximal isotropy subgroups in Figure 12.11 here. We have already explained the 
isotropies 0(2)~ and 0(2) © 7L C 2 of axisymmetric v G spanY^o in (12.121) . Since —id ^ S*0(3) 
acts as (— l) f id on Ve, any element of Vt possesses isotropy at least 7L C 2 = {±id}, if £ is even. 
Therefore any closed isotropy subgroup r„ of T = 0(3) takes the form r„ = K © = 
K U (-K) with K < 50(3). The subgroups of 50(3) are 5*0(2), 0(2), T, O, /, the discrete 
cyclic subgroups Z n of 50(2), and the dihedral subgroups D n of 50(2). Here D n is the 
subgroup of planar rotations and reflections in 0(2) < 50(3) which fix a regular n-gon. 
Similarly T,0,I < 50(3) fix the regular tetrahedron, octahedron, and icosahedron and are 
isomorphic to the alternating symmetric group A±, to the symmetric group 54, and to the 
alternating symmetric group A$, respectively. 

The "octahedral" group 0^ < 0(3) is the group of orthogonal rotations and reflections 
which fix a regular tetrahedron: this time the subgroup of rotations is the tetrahedral group 
T = 0~ H 50(3) of index 2 in 0~. Similarly the dihedral group D§ < 0(3), isomorphic to 
D Q < 0(2), has index 2 subgroup D d % n 50(3) = D 3 < 0(2). 
We summarize the results of this section, following [231 12H E31 E2] . 
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Figure 2.1: TTie lattices of conjugacy classes of isotropy subgroups of the natural irreducible 
representation of 0(3) on £/ie spaces Ve of spherical harmonics, £ = 0, . . . , 4. See 07] / 



Lemma 2.1 (Lauterbach et al.JM ES EfJ 

Consider the irreducible representations ( 12. 8ft of T = 0(3) on £/ie spherical harmonics 
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subspaces Vg; see (I2.10p . 

Then the lattices of conjugacy classes of isotropy subgroups T v on Vg and the real dimensions 
of the fix spaces Fix(T v ) = V?" are given in Figure 12.11 The one- dimensional fix spaces of 
maximal isotropy subgroups are given in Table l2J\ for axisymmetric solutions and all I, and 
for the remaining cases in £ = 3,4. 





0(2) ©Z^ 


0(2)- 




o- 


0®Z C 2 


£>0 


even 


odd 


3 


3 


4 


e G Fix(r„) 


Ygo 


Ygo 




ReF 32 


F 40 + v/5/14r 42 



Table 2.1: Explicit generators e of the one-dimensional fix spaces Fix(T v ) of isotropy 
subgroups T v in certain spherical harmonics representations Vg of T = 0(3). 

The graphs over S 2 of the generators e f the one-dimensional fix spaces Fix(r^) of isotropy 
subgroups r„ listed in Table 12.11 are represented in Figure 12.21 



3 Equivariant branching 

The equivariant branching Lemma 13.11 below was developed by Vanderbauwhede [38] into a 
simple, but strikingly effective tool of equivariant bifurcation analysis; see also [14] . As we 
will see, it is a straightforward adaptation of the classical bifurcation theorem by Crandall 
and Rabinowitz on bifurcation of zeros of maps [T3] : 
We consider O fc -maps, k > 2 

, . G : RxX^Y 

1 j (X,v) i-> G(X,v) 

with Banach spaces X, Y and a trivial branch 
(3.2) G(A,0) = 
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Figure 2.2: The graphs of the generators e of Table [Ql over the sphere S 2 

of zeros. Moreover we assume T-equivariance, i.e. 

(3.3) G(\,yv) = yG(\,v) 

for all 7 in the group T, all A G K and u G X. In particular the observations of Section [2] 
apply to F = G(X,v); see (12. 3p . For any subgroup K < Y let 

(3.4) G K :RxX K ^Y K 

denote the restriction of G to K x X A , where X , Y x are the i^-fixed closed subspaces of 
X,F introduced in (J2TB) . f!2T7j) . 

Lemma 3.1 (VanderbauwhedeJF^) 

In the setting ( I3.ip ~ (j3.4p assume there exists a subgroup K <T and Xq G K suc/i t/iat i/ie 
linearization L K (X) = D X G K (X, 0) zs Fredholm of index zero. Moreover assume 

(3.5) kerL^(Ao) = span {e} 
rea/ one- dimensional, and the transversality condition 

(3.6) <9 A L K (A )e £ Range L K (X ) 
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holds. 

Then the nontrivial solutions of G \X,v) = near (X , 0) form a bifurcating C k ~ 1 curve 

(3.7) s i-> (X(s),v(s)) Glx X k 
which satisfies 

(3.8) A(0) = A , u(0) = 0, u'(0) = e. 

In particular the bifurcating curve intersects the trivial branch (A,0) only at (Ao,0) ; where 
s = 0. TTie isotropy T v ^ is at least K, along the branch. 

Proof: 

Statements (13 .7p . (13. 8p follow directly from [15]. The intersection claim is a consequence 
of v'(0) = e ^ 0. The isotropy claim follows from v (s) G X^; see the example following 
definition (12. 6p of This proves the lemma. M 

4 Center and strong stable manifolds 

We shall now obtain invariant manifold theorems for Equation (ll.ISp which rewrites as 

(4.1) d t v = (v + l) 2 (Av + Xf(v)) 

at an arbitrary equilibrium v*. We regard this as an abstract dynamical system in the 
Banach space X = C^(S 2 ), (3 e (0, 1). and follow the approach of [25] closely. 
Notation: C s , where 5 G R, 5 = [5] + {5}, [5] G N entire part of 5, {5} G [0, 1) its decimal 
part, denotes the space of [<5]-times differentiable functions whose [5]-th derivative is {5}- 
Holder continuous. 

We take X = C 13 and D = C 2+l3 (S 2 ) as the domain of the right hand side. Let v — v — v*. 
We must rewrite the equation in terms of v* and v. This is not difficult since the right hand 
side is just a polynomial in v, u*, Av, Av*. We shall not carry this out in full detail here, but 
shall simply note that the equation takes the form 

(4.2) v t = Av + Q(v), 

where A is just the linearization of the right hand side at u*, and Q is the nonlinear part, 
which satisfies Q(0) = 0, Q'(0) = 0. The operator A is easily computed to be 

(4.3) A = (1 + v*) 2 ( A + A (1 + - 
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To reduce the technical requirements of the present paper, at this stage, we observe that the 
linearization of (ll.lSp or ( 14. 2 p is Morse equivalent to the semilinear parabolic PDE 



(4.4) d t v = Av + f(v), 

again with f(v) = 2v — v 2 /(l + v) as in ( I1.19p . By Morse equivalent we mean that the 
total algebraic multiplicity i(t>*) of the positive spectrum cr + of the linearization at any 
equilibrium of f !4.4j) coincides with that of the linearization A of fll,18p . 

Lemma 4.1 Equilibria of v* of M.18\) . l[4-4\ ) are Morse equivalent. 

Proof: 

By the L 2 self-adjointness of = A + \f'(v*) we notice that all eigenvalues \i of the 
linearization and of the Linearization A = (1 + v*) 2 L* are real. Indeed, Aw = fiw in 
L 2 (S 2 ) implies 
(4.5) 

H / (1 + v*y 2 ww = fx ((1 + f*)" 1 ^, w) = {L*w, w) = (w, L*w) = (w, L*w) = (w, L*w) G R 
Js 2 

for the standard Hermite L 2 scalar product (•, •) on S 2 . Hence /j £E. Moreover, /x = is 
an eigenvalue of if and only if /i = is an eigenvalue of A with equal multiplicity. The 
homotopy of operators 

(4.6) L(t) = (1 + tv,) 2 L, 

with < t < 1 shows the same property. Therefore the linearizations and A are Morse 
equivalent and the lemma is proved. 

Here it is important that the algebraic and geometric multiplicities coincide. We prove this 
indirectly: suppose the algebraic multiplicity of /i is strictly larger. Then there exists a 
nonzero eigenvector w\ with /awi = Aw\ and W2 such that ui\ = (fi — A)w2- Abbreviating 
m — 1/(1 + v*) 2 > 0, this implies 

(4.7) = (m/i - L*) wi 

(4.8) mw\ = (mfi — L*) w 2 - 

Because and /im are self-adjoint this implies 

(4.9) (wi, mwi) = (wi, (m/i - L*)w 2 ) = ((m/i - L*)u>i, w 2 ) = 0. 

This contradiction to positivity of mw\ proves that the geometric and algebraic multiplicities 
of L coincide. txi 
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The little lemma above enables us to derive Morse indices as well as the multiplicity of 
any purely imaginary spectrum from considerations of only the semilinear parabolic equation 
(I4.4p instead of the quasilinear problem (j!.18p . For semilinear semigroup settings, spectra 
and their associated center, center unstable and strong stable invariant manifolds have been 
studies very extensively in the existing literature. See for instance [20], [37], [5], and the 
many references there. In particular, we will invoke center manifolds in our local bifurcation 
analysis of (jl.lBp in Section 5 and 6, to establish normal hyperbolicity and determine the 
Morse indices of nontrivial equilibria ^ of reduced isotropy, which bifurcate from the 
trivial equilibrium v * = at A = Aj = 1(1 + 1). 

Once these Morse indices of normally hyperbolic equilibria Tv* ^ have been determined, 
however, we need to study the remaining solutions v(t) of the original equation (11.181) which 
converge to t>* for t — > oo. These solutions v(t) form strong stable manifolds of u*, i.e. 
manifolds characterized by an exponential decay with bounded exponentially weighted norm 

(4.10) ||i> — f*|| e := sup \\v(t) — v*\\ e et 

i>0 

for small enough e. See [22], [28] for pertinent results in a quasilinear setting. Next, we 
adapt such a result in the setting of [25] . 

As described previously, we work with equation (14. 2 p which is the result of a self-similar 
rescaling at the ODE blow up rate, and centering on an equilibrium v* of the rescaled 
equation. The ambient Banach space is chosen to be X = C^^S 2 ), (3 G (0, 1). The domain 
D of the linear operator A is defined as X D D := C 2+ "(S 2 ). Our approach consists in 
applying Theorem 9.1.7 of [25], and for that we introduce some basic interpolation spaces 
between D and X, and state what they turn out to be in the case of the equation (14. 2 j) we 
are studying. 

Definition 4.1 For a G (0, 1), we define 

D A (a, oo) = {v G X : t m- \\t 1 ~ a Ae tA v\\ x G L°°(0, 1)} 

and 

IMId. 4 Koo) = \M X +[v\d a {*,oo), 

where \v\d a {o.,oo) denotes the L°° norm in time over [0, 1] of the defining function | i 1-£ M.e tA i; |^ 
Furthermore, we recursively define 

D A {1 + a, oo) = {v G D : Av G D A (a, oo)} , 

with 

IMlD A (l+a,oo) = \\ V \\ + II ^11 D A (a,oo) " 
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Remark 4.1 The interpolation spaces defined in 4-1 ore Banach spaces if the operator A is 
sectorial. In our settings, 



(4.11) X = C?(S 2 ), fie (0,1), 

(4.12) XdD = C 2+ P(S 2 ), 

it turns out that 

(4.13) D A (a, oo) = C 2a+ P(S 2 ). 

This result is a straightforward adaptation of Theorem 3.1.12 in l2h^ . here with space variable 
x G S 2 (or any compact manifold without boundary) instead o/iel™. 

In order to state the strong stable manifold theorem, let us introduce the projection on the 
unstable and strong stable eigenspaces respectively. More precisely: Let a (A) denote the 
spectrum of A and take cr+(A) = {/i G cr(A) : Re/i > 0} and cr-(A) = {fi G cr(A) : Re/i < 
0}. Note that o~+(A) consists of finitely many eigenvalues of finite multiplicities associated 
with finitely many eigenfunctions. In the case of the trivial equilibrium v* = 0, these 
are just the spherical harmonics. Take now P + to be the eigenprojection from X onto the 
associated eigenspace and P_ = (id —P+)- The following strong stable manifold theorem is an 
adaptation of [25], Theorem 9.1.8 with the difference that we allow the linearization A to have 
a zero eigenvalue. To compensate for this, we use the spectral gap: let w_ = sup a-(A) < 0, 
and r\ G (0, —uJ). For any time interval / C R, we define the space of functions decaying at 
least as fast as e~ vt by 

C v (I,D A (a+l,oo)) = {g:t^e*g(t)eC(I,D A (a + l,oo))}, 

||fi , lla ) (/.£> A (cH-i,oo) = sup We^gWdl^Aia+l.oo)- 

tei 

Theorem 4.1 For all a G (0, 1), there exists a constant B>\ > and a Lip schitz function ip, 

(4.14) V : 5 (0, h) C P-(C 2a+ P(S 2 )) -> P + (C 2(Q+1)+/3 (^ 2 )), 

differentiate at with ip'(0) = 0, whose graph is the local strong-stable manifold for Equa- 
tion (14.21) . In other words: 

For all v G Graph(ijj), equation \4-2\ ) has a unique solution v G C v ([0, oo(, C 2( - a+1 ' + ^ (S 2 )) 
such that \\v\\c v ([o,oo),c^+W(s 2 )) < B i- 

And conversely there exists a constant b\ > such that if equation \4-ty has a solution v G 
C v ([0,oo),C 2 ( a+1 ^(S 2 )) with \\v\\ Cn {[v,oo),cWW{s*)) < Bi and \\P_{v{0))\\c^{S*) < h, 
thenv(0) G Graph(ip). 
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This is essentially part (ii) of Theorem 9.1.8 of [25] or part (i) of Theorem 2.4 of [27], i.e. 
the strongly stable manifold parts of those theorems. Those theorems apply to a dynamical 
equation of the form (14. 2p . where A is sectorial and the nonlinearity Q may be regarded 
as a C°° mapping from a neighborhood O C D^(a + 1, oo) into DA_(a, oo) and in addition 
satisfies Q(0) = 0,Q'(0) = 0. Taking the interpolation spaces as generated by C^{S 2 ) and 
C 2+l3 (S 2 ) as described above, these conditions hold in our case. However, the results from 
[25] , [27] both assume that there is no null eigenvalue, which we do not assume. It can 
nonetheless be checked that one retains a strongly stable manifold theorem even when an 
eigenvalue vanishes as long as one utilizes the spectral gap (0,u;_) correctly. 
In order to discuss the heteroclinics between the branches of equilibria (see Section 7, Remark 
17.11) . we shall prove a center-unstable manifold theorem. 

To carry out this discussion we recall that the operator A : D — > X is sectorial. Again, it 
is important that the nonlinearity Q is a C°° mapping Q : C 2{a+1)+/3 (S 2 ) C 2a+f3 (S 2 ), 
which in addition satisfies Q(0) = 0,0/(0) = 0. Let A + = A P+ ^ C 2 a +p^ S 2^ and A_ = 
y4p_( C 2( a +i)+/3( S 2)). With x{t) = P+v(t) and y{t) = P-v(t) we make the following decom- 
position of our equation (11.181) : 



For the center-unstable manifold theorem we shall first make a small modification of the 
system KTty - Km . Let p : P + {C 2a+p (S 2 )) R be a C°° cutoff function such that 

< p(x) < 1, p{x) = 1 if || ^Hc'2a + /3(-5'2') 1/2, p{x) — if ||x|| C 2a + / 9(52) > 1- 

The existence of such a function is no problem since P + (C 2a+l3 (S 2 )) is finite dimensional. 
For q > small we consider instead the system modified by cutoff 



where Q + (x,y) = P + Q(p(x/q) + y) and Q-(x,y) = P-Q(p(x/q)x + y), and q is called the 
cut-off parameter. 

note that the above construction can be carried out in such a way that it is equivariant 
with respect to 0(3) in the case that the original equilibrium is trivial i>* = 1. Indeed, it 
is certainly the case that our original equation, as well as the spaces X + = P + (C 2, ^(S' 2 )), 
X_ = P_(C 2, ^(S' 2 )), are equivariant under this group action. Hence we must only check that 
we can arrange for the modified equations (I4.17p - (l4.18p to be equivariant as well. In order to 
do this, we use an isometry (using the L 2 norm on X + ) F : X + :— > W n to identify X + with 



(4.15) 
(4.16) 



x'(t) = A+x(t) + P+Q(x(t) + y(t)) 
y'(t) = A-y(t) + P-Q(x(t)+y(t)) 



(4.17) 
(4.18) 



x'(t) =A + x(t)+Q + (x(t),y(t)) 
y'(t)=A_y(t) + Q.(x(t),y(t)) 
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M. m . Thus, arranging that the cutoff function p is such that p o F is spherically symmetric, 
we can ensure that the modified equations (14.1 7p - (14.181) are 0(3) equivariant. 

Theorem 4.2 Let now A and Q be as in the first paragraph of this section and consider 
equation fl^.i^ ) as well as the modified equations fl^.i7| )- P.i&| ) constructed with a cutoff p as 



in the remarks of the previous paragraph. 

In particular A is sectorial and admits finitely many unstable eigenvalues of finite multi- 
plicities. The nonlinearity Q is assumed to be a mapping C QC (0,C 2a+l3 (S 2 )), where O is a 
neighborhood of the origin in C 2( - a+1 ^ + P (S 2 ). 

Then there exists a cutoff parameter qi such that if q < g 1; then there is a continuously 
differentiate function E : P + (C 2a+l3 (S 2 )) — > P_C 2( - a+1 ^ +l3 (S 2 )) whose graph is invariant for 
the system l\4-17\ )- [4.18\ ). GraphE is called a local center-unstable manifold. In addition 

E'(x)(A + x + Q + (x,E(x))) = A.E(x) + Q-(x,E(x)),x e P+{X). 

In the case of the trivial equilibrium, E is equivariant with respect to 0(3), and so is its 
graph. 

The function E can be made C k by choosing a possibly smaller bound q^ for the cutoff 
parameter. 

This theorem can be found in [25]. For a similar theorem see the nice succinct work of 
Mielke gB]. 

Of course the finite dimensional center-unstable manifold contains a center manifold which is 
obtained by using the projectors P an d Ph onto the center eigenspace associated to the zero 
eigenvalue and the "hyperbolic eigenspace" of infinite dimension associated to the nonzero 
eigenvalues respectively. With x{t) = Pov(t) and y{t) = Phv{t) we make the following 
decomposition of our equation fll . 18j) : 

(4.19) x'(t) = P Q(x(t)+y(t)) 

(4.20) y'(t) = A h y(t) + P h Q(x(t) + y(t)) 

This proves the following theorem. 

Theorem 4.3 Let O be a neighborhood of the origin in (7 2 ( a + 1 )+/ 3 (S 2 ) in which the non- 
linearity Q is continuously differentiate, Q E C x {0, 2< - Q+1 ' )+/3 (S' 2 )). The linear part A is 
sectorial and admits only a finite number of nonnegative eigenvalues of finite multiplicity, 
Q(0) = Q'(0) = 0. 

Then there exists a cut-off parameter q\ such that for < q < q\, there exists a Lipschitz 
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continuous function $ : Po(C 2a (S 2 )) — >• P/ l (C 2 ^ a+1 ^(5' 2 )) whose Graph is invariant under the 
flow of the system modified by cut-off 

(4.21) x'(t)=Q (x(t),y(t)) 

(4.22) y'(t)=A h y(t) + Qh(x(t),y(t)), 

where Q (x,y) = P Q(p(x/q) + y) and Q h {x,y) = P h Q(p(x/q)x + y). 
Graphs is called a (local) center manifold. 

5 Anisotropic self-similar blow-up 

In this section we apply the equivariant branching Lemma 13.11 to obtain anisotropic self- 
similar blow-up solutions 

(5.1) u(r,p)= 2 (v + 1) 



2 \r 

with isotropically prescribed scalar curvature R{r) = (A + 2)/r 2 , as we announced in the 
introduction. By (11. 1 8 j) this requires v = v(p) to be a nonhomogeneous solution of the 
equilibrium equation 

(5.2) = Av + Xf{v) 

for the blow-up shape v(p), p G S 2 . We also recall the nonlinearity 

(5.3) f(v)v - l -v 2 /{v + 1) 
and the coefficient 

(5.4) A = r 2 R(r) - 2 

Anisotropic, i.e. nonhomogeneous, solutions v of (I5.2p will arise by 0(3) symmetry breaking 
bifurcations with maximal isotropy and one-dimensional fix spaces, as indicated in Lemma 

ED 

To apply the equivariant branching Lemma 13.11 we first rewrite (15. 2 j) as 

(5.5) = G{\, v) := v + (A — l)-\v + \f(v)). 
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For Banach spaces X = Y we could choose the Hilbert space setting X = L 2 (S 2 ). For 
technical consistency with our dynamics approach of Holder type, above, we work in 

(5.6) X = Y := C?(S 2 ) 
instead. 

We notice immediately, that G is a compact perturbation of identity. Indeed the solution 
(A — l)~ 1 g of the linear Poisson equation Av — v = g on S 2 is a bounded linear operator 

(5.7) (A- l)- 1 : C P {S 2 ) ^C?(S 2 ), 

and the embedding C 2+l3 (S 2 ) C^{S 2 ) is compact by Arzela-Ascoli. 

Obviously G is C k , for 1 + v > strictly positive and any k > 0, by analyticity of / on 
this subdomain of R, and of C Q,a . A trivial branch G(X, 0) = of homogeneous, 0(3)-fixed 
solutions arises from /(0) = 0. Orthogonal invariance of the Dirichlet integral | J |Vt> \ 2 dp 
under the representation (Q( r y)v)(p) := v( r y~ 1 p) of f)2.8p for 7 G T := 0(3) implies 0(3)- 
equivariance f)2.3p of (A — 1) _1 on v G C ' a (S 2 ). Obviously 0(7) also commutes with the 
point evaluation v H- f(v) in /3 (S' 2 ). This proves T-equivariance of G and establishes 
assumptions fl3.ip - fl3.3l) of the equivariant branching Lemma [3.11 
Consider the 0(3)-equivariant linearization 

(5.8) L(X)v := D v G(\,0)v = v + (A-iy\v + \f'{0)v) = 

= u + (l + A)(A-l)- 1 v 

at the trivial branch G(X, 0) = next. Again L(X) is a compact linear perturbation of 
identity on X = Y = C^(S 2 ). Hence L(X) is Fredholm of index zero on X, and on any 
closed linear and L(A)-invariant subspace of X. In particular L K (X), the restriction of L(X) 
to the subspace X K = Fix(K) of i^-fixed elements in X, is Fredholm for any subgroup K 
of T = 0(3). 

By standard spectral theory of the Laplace-Beltrami operator A on the 2-sphere S* 2 we obtain 
nontrivial kernel 



(5.9) {0} ^ kerL(A) = V e 
if and only if, 

(5.10) \ = \ e ■=£(£+ l) 

for some £ = 0, 1, . . .. The resulting kernels are the L 2 -complete spherical harmonics sub- 



spaces V e oiX = C°' a (S 2 ); see fl22}, (12TTUD in Section H 
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To check the transversality condition (13. 6p of the equivariant branching Lemma 13.11 we cal- 
culate 



(5.11) D x L(\)e = (A - l) _1 e 



+ 1) 



e. 



By L 2 self-adjointness of A, however, e e kerL(A) is L 2 orthogonal to range L(X) and, a 
fortiori, orthogonal to any subspace range L K (X). This proves transversality condition ( 13. 6p . 

Theorem 5.1 

Let £ = 1, 2, 3, . . . and let K < 0(3) be an isotropy subgroup of 0(3) such that 

(5.12) V e K = Fix(K) = span {e} 

is one-dimensional, for the representation (g( , y)v)(p) := v( r y~ l p) of 0(3) in the spherical 
harmonics space Vi of dimension 2£ + 1; see A2.8\) , 112. 10\) . 
Then there exists a unique smooth local curve 

(5.13) s ^ (X(s),v(s)) G M x C\S 2 ) 

of solutions (A, v) of the self-similar blow-up equation A5.2\) such that 

(5.14) A(0) = X e = £(£ + 1), u(0) = 0, v'(0) = e 
and the isotropy ofv(s) is exactly K < 0(3), for s ^ 0: 

(5.15) T V{S) = K^0(3). 

In particular v(s) is nonhomogeneous and hence anisotropic, for s / 0. 
Proof: 

By the arguments above and assumption ( 15. 15ft . the equivariant branching Lemma [37X1 asserts 
the existence of a local branch (X(s),v(s)) elx X K , of class C k for any k > 1. Smoothness 
follows, for the k = 1 branch and any s ^ 0, by direct implicit function theorem. At s = 
smoothness follows by uniqueness of the branches, even if the s-neighborhoods had shrunk 
to zero for k — > +oo. 

To check for the precise isotropy T v ^ = K is a little more subtle. In Section |5] we have 
already remarked that isotropies K with one-dimensional fix spaces are maximal in their 
irreducible representation. This proves 

(5.16) r e = K ^ 0(3) 
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for the eigenvector e which spans V e . 
To extend this observation from e to 

(5.17) v(s) = se + o(s), 

for small enough s, we argue indirectly: suppose there exists a sequence 7^ s, — > such 
that 



(5.18) #<r„ (s ,) <r = o(3). 

In r = 0(3) this leaves only finitely many choices for r\,( Si ). Indeed the only infinite ascending 
series of closed subgroups are generated by the rotations over 2n/k in 50(2) < 50(3). 
This case is easily excluded, by direct inspection of Fourier series in the azimuthal latitude 
(f. Without loss of generality we may therefore assume 

(5.19) K < T v{Si) = ^ < T = 0(3) 
does not depend on i. 

Since v'(0) = e we have v(s.i)/si — > e for 7^ Sj — > 0. Continuity of the representation g 
therefore implies 



(5.20) T v{si) =T 00 <T e = K 

which contradicts our indirect assumption (15.181) . This argument, which is a special variant 
of the virtual isotropy proposition in [16] , completes the proof of (I5.15P and of the theorem. 

cxi 

We now apply Theorem 15.11 on 0(3) symmetry breaking of the trivial equilibrium v = 
in (15.51) to the maximal isotropy subgroups K of the spherical harmonics representations 
of Lemma [2.11 We address the axisymmetric cases K > 50(2), for all £, and the discrete 
isotropics K, for £ = 3,4, separately, in Corollaries 15.11 and 15.21 below. 

Corollary 5.1 

Let i > and consider axisymmetric fix-groups 



(5.21) K :-- 




for odd i 
for even 



Then the conclusions of Theorem \5. 1\ hold and provide a unique smooth local branch (A(s), v (s)) 
of nontrivial solutions bifurcating from (A(0), v(0)) = (Xe, 0) at Xi = £(£ + 1). The solutions 
v(s) at X = X(s) possess isotropy r v r s -\ = K for s 7^ 0. In particular they are nonhomogeneous 
but axisymmetric, i.e. independent of the azimuthal angle ip on S 2 . 



24 



Proof: 

Lemma 12.11 implies that the one remaining assumption ( I5.12p of Theorem 15.11 holds for K 
from fl5T2T|) : 

(5.22) dim Fix (K) = dimV e K = 1. 

This proves the corollary. 

IX] 



Corollary 5.2 

Let £ = 3 and K e {0~, Df], or let i = 4 and K = O © Z c 2 . 

Then the conclusions of Theorem \5. 1\ hold and provide a unique smooth local branch (A(s), v (s)) 
of nontrivial solutions bifurcating from (A(0), v (0)) = (A^, 0) at Xi = £(£ + 1). The solutions 
for v(s) possess isotropy r„( s ) = K , for s ^ 0. In particular they are nonhomogeneous with 
hexagonal or octahedral symmetry as given by K . 

Proof: 

The proof of the previous corollary applies verbatim. 

IX] 

In conclusion we obtain self-similar but anisotropic blow-up in the isotropic scalar curvature 
equation, with remaining octahedral, hexagonal, or axisymmetry as specified by the isotropy 
K above. 

Of course our result applies equally well to all maximal isotropy subgroups K of 0(3) on 
Ve, as long as the spaces Fix(K) of infixed vectors in Vi remain real one-dimensional. See 
[T3l [23] for lists of these cases. 



6 Anisotropic asymptotically self- similar blow-up 

In this section we address dynamic aspects of the self- similarly rescaled scalar curvature 
blow-up problem 

(6.1) (l + v)- 2 d t v = Av + Xf(v); 

see (I1.18p . More specifically we study the Morse indices of the 0(3) symmetry breaking 
bifurcations (X(s),v(s)) from the trivial homogeneous solution branch (A,t>) = (A,0) at the 
bifurcation points (A(0),c(0)) = (A^,0) with 

(6.2) X e = 1(1 - 1) 



25 



for I = 1, 2, 3, . . .. The nonlinearity f(v) takes the specific form 

(6.3) f(v) = v - ^7(1 + v) = v - l -v 2 + ^ 3 T ■ ■ ■ 
In particular the first Taylor coefficients are 

(6.4) /(0) = 0, f(0) = 1, f (0) = -1, f(0) = 3. 

We give a complete account of the low-dimensional bifurcation diagrams i = 1,2,3,4 in 
Theorem 16.11 based on Corollary 15.11 and 15 .21 We omit the exceptionally trivial case I = 
0, \i = 0, which is linear and only features vertical bifurcation of spatially homogeneous, 
isotropic equilibria v = const. 

Our local analysis starts from the standard equivariant center manifold theorem with equiv- 
ariance group T = 0(3) for the Morse equivalent but semilinear parabolic PDE 

(6.5) d t v = Av + Xf(v). 

Adding the equation A = 0, artificially, we obtain a trivially extended semiflow S x (t) of 
class C k on the extended phase space Rx(Ifl{l+t)> 0}). At bifurcation points (A, 0) = 
(A^, 0), A^ = £(£ + 1), we immediately obtain a reduced flow of (\6.5\i on the center manifold 
tangent to Vi in parametrized and equivariant form: 

(6.6) v c = $(\,v c ), v c eV e . 

In Lemma 16.11 we first compare this flow, truncated to order 2 or 3, with the detailed 
discussion of the generic 0(3)-equivariant bifurcations on Vi in [231 1X3] . Specifically, we 
establish nondegenerate transcriticality and the sub- or supercritical nondegenerate pitchfork 
type of the maximal isotropy branches of Corollaries 15. 1[ I5.2[ respectively. We call the 
bifurcation of (\(s),v(s)) at s — 0, (A(0),w(0)) = (A£,0) nondegenerate and 

trans critical, if A'(0) ^ 

(6.7) 

pitchfork, if A'(0) = ^ A"(0). 

Because these branches of equilibria must reside in the local center manifold of f)6.5p . they 
establish nondegeneracy of the respective quadratic and cubic equivariants of x c i— > Q(\e, x c ), 
and thus determine the resulting unstable dimensions of the bifurcating equilibria v(s) inside 
the center manifold. This little trick avoids very messy direct computations of, and inside, 
the center manifolds. 

In a second step, in Proposition 16.21 we then establish normal hyperbolicity of the associ- 
ated group orbits Tv(s) = T/T v ^ of equilibria v*. We determine the codimensions of the 
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local strong stable manifolds W ss = graph\l/ ss , for each of the symmetry breaking branches 
specified in Theorems 15. 1\ 16.11 and Corollaries 15. 1\ 15.21 Since 



(6.8) 



v(-,t) -> 



converges exponentially in C 3 , on graph \I/ SS , the solutions in graph ty ss provide anisotropic 
and only asymptotically self-similar solutions u with blow-up at r = 1 of the isotropic scalar 
curvature R{r) = (A + 2)/r 2 . 

Lemma 6.1 

Consider the symmetry breaking branches (Xk{s),vk{s)) of equilibria of Equation (16. ip . [6751 
with isotropy K which bifurcate from (\k(0),vk{0)) = (Xe, 0), A^ = £(£ + 1), £ = 1, 2, 3, ... , 
as established in Corollaries 15.11 and 15.21 

Taen i/ie bifurcations of the axisymmetric solutions K > 5*0(2) of Corollary 15.11 are non- 
degenerate trans critical, for all even £ > 0. TTie remaining solution branches at £ = 3,4 
established in Corollary 15.21 are 



At £ = 3,4 the pitchfork curvatures A'^(O) o/ £/ie branches with maximal isotropy K are 
given in Table lfTM Here we use v' K (0) = e# to normalize the parametrization of all branches 
(\k(s),Vk(s)) by s such that 



holds for the L 2 scalar product (•, •) on S 2 with the normalized eigenvector of the space 
V e K of K- fixed vectors in V^. 

Proof: 

We first calculate the derivatives of A(s) at s = 0, in the general setting of the symmetry 
breaking bifurcation Theorem 15.11 We then specify to the axisymmetric cases K > SO (2), 
and to the maximal isotropics K in the Ve representations of T = 0(3). 
We recall the general setting 



(6.9) 



nondegenerate transcritical, for £ even 
nondegenerate pitchfork, for £ odd 



(6.10) 



(e K ,v K (s)) 



= s 



(6.11) 



G(X, v) 
L(X)v = d v G(X,0)v 

f(v) 



v + (A-ir 1 (v + Xf(v)) = 

t; + (A-l)- 1 (l + A)t; 

1 1 

v - -v 2 (v + ly 1 = v - -v 2 + 
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isotropy K 


0(2) 


0(2) © Z£ 





© z c 2 


^6 


£ = 1 


e K 


2 10 










A (UJ 


3 

5tt 










£ = 2 






- 1 20 








A'(0) 












£ = 3 


e*- 


^30 




\/2Rer 3 2 




\/2ReF 3 3 


A"(0) 


2954 

7157T 




1050 

1437T 




665 
286tt 


£ = A 






^40 








A'(0) 




243 
2OO20F 




9^ 
2860F 





Table 6.2: Exact values of normalized eigenvectors spanning isotropy fix spaces V e and 
of first nonvanishing derivatives \' K (0), A'^-(O) for symmetry breaking equilibrium branches 
(A^(s), vk{s)) of (16. II) . (16.51) . with maximal isotropy K. 
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from (JE3D, A53J) . ( E5D - By definition 

(6.12) Fix (A) n ker L(X e ) = Fix (A') nV t = span {e^}, 

where we choose to be an L 2 -normalized unit vector. See Theorem 15.11 and f)5. 12[) — (1 5. 14p . 
where dependence on the chosen isotropy A was suppressed. For the explicit entries in 
Table 16.21 see [T5] and [T7]. The case e K = Y i0 for A > SO (2) is obvious: the fix space 
generator ex £ Vi is independent of the azimuthal angle </? and hence normalized to ±Y£o($) 
in the spherical harmonics space Vg. 

From Section[2]we recall that —id G Z2 acts on by multiplication with (— 1) £ . In particular 
—id ^ A for any nontrivial isotropy A on V% and odd £. Conjugation by —id G Zj, on the 
other hand, leaves K invariant. Uniqueness of the bifurcation branch (\k(s),v(s)) and the 
normalization (I6.10p therefore imply 

(6.13) Xk{s) = \ K (s), v k (-s) = -v K (s) 
for any isotropy K and odd I. In particular 

(6.14) A' A -(0) = 

and all such symmetry breaking bifurcations are pitchfork candidates. 

To calculate the first derivative of A(s) at s = 0, in general, we suppress K again and 

differentiate G(X(s),v(s)) = twice with respect to s: 

(6.15) = d x GX' + d v Gv' 

(6.16) = d x G\" + dlG{\') 2 + 2d x d v GX'v' + d 2 v G{v'f + d v Gv" . 

At s — 0, v(0) = all A-derivatives of G vanish because G(X, 0) = along the trivial branch. 
Self-adjointness of A implies orthogonality 

(6.17) (e,d v G) = 

in L 2 (S 2 ) at s = where d v G = L(X(). Taking the scalar product of (16.141) with e at s = 
implies 

(6.18) = 2A'(0)(e, (A - l^f^e) + (e, (A - I)" 1 /" (0)e 2 >. 

With Ae = — A^e, /'(0) = 1, f"(0) = —1 and the L 2 -normalization (e, e) = 1 we obtain 

(6.19) A'(0) = ^(e,e 2 >. 
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We claim next that 



(6.20) A^(0) ± 

is nondegenerate transcritical, for all axisymmetric solutions K = 0(2) © 7L C 2 and even 
2, 4, 6, We use [I], 27.9.1 to explicitly calculate (e, e 2 } = j (Y eo ) 3 as follows 

S 2 



(e, e 2 ) 



(21 



(6.21) 



47T 



\0 0, 

?f 4- 1 

(eeo o\££io) 2 




We have used I even in the last equality; for odd I we encounter a rather circuitous con- 
firmation of A'(0) = 0. The alternating sum from /c = 0to/c = £/2is nonzero because 
the binomial coefficients increase strictly for k — 0, . . . ,£/2. This proves nondegenerate 
transcriticality (I6.20p for axisymmetric solutions and even i. For example we obtain the 
entry for i = 4, K = 0(2) © %\ in Table | 



We consider the octahedral case, £ = 4, K = 0©Z§, e = e K = \Y 40 +2y/5/14 Re F 44 ^ 

next. Up to positive coefficients we obtain 

(6.22, <*«?> = / (I^(no + ^(y„ + ^))) 3 ^^^0. 

s 2 

The last equality can be verified, e.g., by Mathematica or similar programs. Alternatively it 
follows by hand, again, using [TJ, 27.9.1. 

It remains to determine the three pitchfork curvatures A'^(0) for £ = 3, K G {0(2)~, 0~, D$}, 
e = e#. A similarly tedious but crucial calculation was already performed in [17], based 
on the generators e = of the isotropy fix spaces from [23| [13]. We differentiate the 
second derivative (I6.16p . once again, with respect to s at s = 0, using A'(0) = and <9f G=0 
for all k: 

(6.23) = 3d x d v GX"v' + d 3 G(v') 3 + 3d*Gv'v" + d v Gv"'. 
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We test against the L 2 unit vector e = v'(0) as before, and insert the derivatives of G, f to 
obtain 

(6.24) = 3A" + 3X e (e, e 3 ) - 3A £ (e 2 , v") 

at s = 0. Evaluating (I6.16p . once again with d\G — 0, A' — 0, v' — e, at s — we observe 

(6.25) (A + Xt)v"(0) = +\ e e 2 . 

Because v"(0) e X K is also L 2 -orthogonal to the one-dimensional kernel V e K = span {e} of 
(A + Xi) in X , by normalization f)6.10p . we can write w"(0) uniquely as a finite sum 

(6.26) v"(0) = V — c e , m , Y e , m , 

\m'\<e' 

where cg/ m i are the coefficients of the spherical harmonics expansion 

(6.27) e 2 = ^2 c t'm> Yi'm'- 

\m'\<l' 

Summarizing (16.24j) - (16.27j) we have the explicit expression 



(6.28) A"(0) = A, Y, 



Af — Af 
\m'\<t> 



I 1 2 

The coefficients c# m / have been tabulated in [17] and are consistent with the detailed 
expansions 

2 in.. i 



(6.29) ei(2) _ = -^Y M + - Yx 
for I — 1, and 

e l{2)- = 2^ ( Y oo + ^ Y 20 + n y 4o + *eo y 

(6.30) = ^(y 00 -iy 40 +^R e y 44 + II lo=y 60 + l^Reya4) 

e l% = 27^ (- y oo + it ^20 - n y 4o + ^Tjj ^60 + ^= Re F 6( 

for £ = 3. Insertion into f!6.28j) proves the remaining i = 3 entries of Table [6721 and completes 
the proof of the lemma. 
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CXI 

The bifurcating branches (Xk{s),vk{s)) of Lemma 16.11 are contained in any local center 
manifold graph f c at A = A^, v — as equilibria of the reduced ODE v c = $(A,t> c ); see 
(16.61) . By T-equivariance, the bifurcating equilibria Ur-(s) of isotropy K, for s ^ 0, generate 
group orbits Tvk(s) = T/K of equilibria which, likewise, appear as C fc manifolds in Graph 
\l/c. The linearized semiflow 

(6.31) d t w = Aw + \f'(v K (s))w 

of (16. 5p at t>x(s) therefore possesses a trivial eigenvalue of geometric multiplicity dimT/K. 
This trivial part of the spectrum is likewise contained in the center part a c which arises by 
perturbation of the eigenvalue of algebraic and geometric multiplicity dim Vi = 2£ + 1 of 
the linearized semiflow (16.311) at A = A^, v — for s — 0. By standard linear perturbation 
theory this center part a c (B x (s)) can be calculated from the linearization of the reduced 
ODE v c = $(A,v c ). 

The remaining non-center part cr^ of the linearization 16.311 is strictly hyperbolic, and is 
inherited from the strictly hyperbolic part at s = 0. This decomposition leads to the following 
proposition. 

Proposition 6.2 

Let (Xk(s),Vk(s)) bifurcate from A = A^, v = as in Lemma \6J\ above and consider small 
enough \s\ ^ 0. 

Then the manifold Yvk{s) of equilibria is normally hyperbolic for the semilinear semiflow 
A6'.5\) if, and only if, it is normally hyperbolic in the center manifold. The strong unstable 
dimension i = i(vk{s)) in X = C^{S 2 ), alias the Morse index, then relates to the strong 
unstable dimension i c within the center manifold as 

(6.32) i = f + i c . 

Normal hyperbolicity, with the same Morse index i(vk(s)) as given by A6.3fy) , likewise holds 
true for the original quasilinear semiflow U.1S\ ), ), Ii6.1\ ) at the same equilibria. 

Proof: Morse equivalence lemma [6J] reduces the quasilinear case to the semilinear case. 
By standard center manifold reduction of the semilinear case, it only remains to prove that £ 2 
is the dimension of the strong unstable part of the linearization 16 . 3 ll at s = 0, A = A^, v — 0. 
Since the strong unstable eigenspace of Re spec< is Vq © ... © Vi-i with dim Vp — 2£' + 1, 
this latter claim is obvious since antiquity, and the proposition is proved. 

ex 

We remark that the strong unstable dimension i in (16.321) is also the codimension of the 
center-stable manifold of (\k{s),vk{s)). It coincides with the codimension of the stable 
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set of the group orbit Yvk(s), which is foliated by the strong stable manifolds over the 
base manifold Y/Yk of equilibria. Although these claims which are standard for semilinear 
semiflows, also hold in the quasilinear semigroup setting, we do not pursue such details here. 

Theorem 6.1 

Consider the branches (\k{s),vk(s)) of nontrivial equilibrium solutions of U.18\) which bi- 
furcate at A_r-(0) = A^, vk{0) = 0, £ = 1,2, 3, 4, and possess maximal isotropy K . 
Then the bifurcating group orbits Yvk{s) = Y/K of equilibria are normally hyperbolic. Their 
strong unstable dimensions are listed in Table \6.3[ . The local bifurcation diagrams are sketched 
in Figure \6.1\ 

£ = £ = 1 £ = 2 £ = 3 £ = 4 




Figure 6.1: local transcritical and pitchfork bifurcation branches (Xk{s), vk{s)) from Afc(0) = 
Xe = £(£ + 1), vk{0) = with maximal isotropics K. Numbers in brackets indicate codimen- 
sions of their normally hyperbolic group orbits. 

Proof: 

As we have seen in Proposition 16. 2\ normal hyperbolicity can be checked within the center 
manifold. For £ < 4, as we consider here, the local stability analysis of bifurcations with 
maximal isotropy has been performed in [231 US] i n terms of lowest order equivariant poly- 
nomials of quadratic and cubic order, respectively, for even and odd £, and under certain 
not overly explicit assumptions on genericity. Closer inspection shows that genericity, in our 
cases, amounts to the bifurcation slopes X' K {0) and curvatures A'^-(O) to be nonzero and, for 
£ = 3, to avoid the absence of the nontrivial cubic invariant. 
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Nonzero slopes and curvatures have been proved in Proposition 16.21 To address cubic non- 
degeneracy for I = 3, we use [T5] , Table 6.2 from which the bifurcation curvatures derive 
as 

{2(18a-/3) If = 0(2)~ 

8(45a-2/3) for K = Df 
16(10a-/3) # = 0". 

Here the general cubic equivariant is written as ac(t> c ) + /3|u c ||u c in [13], with a suitable 
equivariant vector polynomial c(v c ) of the center manifold variable v c G = R 7 . 
The three bifurcation curvatures A^-(O) of Lemma [6.11 and Table consistently overdeter- 
mine the two coefficients a, (3 as 

287 „ 812 

(6.34) a = — , (3 



1430vr ' 143vr 

In particular a ^ 0. By [13], Tables 6.2, 6.3, this establishes normal hyperbolicity of all 
bifurcating branches (Xk(s), Vk (s)), for small |s| 7^ 0, and the strong unstable dimensions 
listed under i c in Table 16.31 These results persist under perturbations by higher order terms; 
see H31, Theorem 5.2. 



For illustration consider the axisymmetric case K = 0(2) _ for I = 3. By Lemma [6. 11 Table 
16.21 we encounter a subcritical pitchfork, A^-(O) < 0, which contributes a simple unstable 
eigenvalue ix$ > to the linearization at (\k(s), i>k(s)) in the center manifold. By [13], 
Table 6.2, the four remaining directions transverse to the two-dimensional group orbit 

(6.35) Tv{s) = r/T v{s) = 0(3)/0(2)- = SO(3)/SO(2) = S 2 

in Vg = MJ produce two further unstable eigenvalues ^1,112 > 0, for cubic coefficients 
This proves i c = 3 and, by Proposition 16.21 

(6.36) i = i c + £ 2 = l2 
as claimed in Table 16. 3[ 

The remaining entries in Table 16.31 are obtained similarly, with group orbits 

Tv = 0(3)/ (0(2) © Z c 2 ) = RP 2 

in the axis symmetric cases £ = 2,4 and with 3-dimensional homogeneous spaces Tv for the 
discrete subgroups K. 
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The local bifurcation diagrams of Figure 16.11 combine the information of Table 16.21 on the 
bifurcation types and directions with the unstable dimensions % of Table 16.31 This proves 
the theorem. 

cxi 

We can now combine the results of the previous chapters to address normal hyperbolicity 
and the strong stable manifolds of the symmetry breaking equilibrium branches (X(s),vk(s)) 
of 11.181 which arise in Lemma 16. 1\ Theorem 16.11 and Table 16.21 



£ 


isotropy K 


type 


direction 


i c 


i 


1 


0(2)- 


pitchfork 


sub 


1 


2 








sub 


2 


6 


2 


0(2) ©Z^ 


transcritical 


super 


1 


5 




0(2)- 


pitchfork 


sub 


3 


12 


3 


o- 


pitchfork 


sub 


1 


10 




u % 


pitchfork 


sub 


4 


13 








sub 


4 


20 


4 


0{2)@Z 2 C 


transcritical 


super 


4 


20 








sub 


5 


21 




0~®Z 2 C 


transcritical 


super 


2 


18 



Table 6.3: For bifurcations (X(s),v K (s)) from \ K (0) = \ e = £{£ + 1), v K (0) = with 
maximal isotropies K and normally hyperbolic group orbits, we indicate the bifurcation 
type, direction, strong unstable dimension i c within the semilinear center manifold and 
strong unstable quasilinear dimension i in X = C^(S 2 ), alias the codimensions of the strong 
stable manifolds of the group orbits Tw^O) for s ^ 0. 



7 Conclusions and outlook 

The existence of a manifold Tv k = T/K of nontrivial equilibria of the self-similarly rescaled 
equation fll,18p . as well as their strong stable manifold W ss (Tvk) has been studied above 
in the context of equivariant bifurcation theory. The normal hyperbolicity of the manifolds 
Tvk has been established in Theorem 16.11 By Theorem 15. 1\ we conclude the existence 
of solutions v of the self-similarly rescaled equation fll . 18j) which converge to a branching 
equilibria v (t, ■) — > vk as t — > +oo. In the original equation II .41 this behavior corresponds to 
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an asymptotically self-similar blow up as r — > r\ = 1 with 

(7.37) u(r,p)=(^(l/r-l)^ * (v{-\ log(l - r),p) + lj , p G S 2 , 

with isotropic scalar curvature R(r) = 
Remark 7.1 (heteroclinic orbits) 

Equivariant bifurcation theory does not only provide a bifurcation diagram of branches of 
equilibria of the self- similarly rescaled equation $1.18\) as given in Figure \6.1\ As Theorem 
provides a local center manifold, we deduce the existence of heteroclinic orbits connecting 
those equilbria. Following Lauterbach et al. |H1 [73] /. recalling the fact that our equation 
satisfies the genericity condition required there as proved in the previous section (see proof 
of Theorem (OOP ), and given the signs of the first nonvanishing derivatives in Table \ 6.^ 
we conclude the existence of heteroclinic connections between the following branches of the 
bifurcation diagram \6.1\ as summarized in the following Table \7-4\ where arrows — > denote 
the existence of heteroclinics orbits between the branches of equilibria. 



t 


type 


A < A^ 


A > A^ 


1 


pitchfork 


0(2) -> 




2 


transcritical 


0(2) ®l? c ^0 


0^O(2)©Z2 


3 


pitchfork 


D% 0(2)- 0- -> 




4 


transcritical 


0(2) © Z c 2 <- © Z 2 C -> 


© Z 2 <r- 0(2) © Z 2 



Table 7.4: Heteroclinic orbits between bifurcating branches. 

_ i 

The self-similar rescaling by the ODE blow up rate (1/r — 1) 2 can be also seen as a Poincare 
"compactification" that allows us to describe the blow up behavior as the convergence to 
an equilibrium in the sphere at infinity. To be more precise, the Poincare compactification 
projects centrally the ambient Hilbert space L 2 (S 2 ) into the unit hemisphere % of L 2 (S 2 ) xl 
(which is a ball of the same dimension as L 2 (S 2 ) itself), in such a way that infinity is projected 
on its boundary, the equator S, alias the sphere at infinity, whose elements will be denoted 
by x ^ L 2 (S 2 ), Hxll = 1 in the following. See Figure 17.21 for a sketch of the Poincare 
compactification and [19] for details. 

This projection is called compactification for historical reasons - it has been developed by 
Poincare for ODE's - but in the case of function spaces such as L 2 (S 2 ), the resulting Poincare 



3(3 




£ 



Figure 7.2: Poincare compactification: each element u G L 2 (S 2 ) is identified with (u, 1) G 
L 2 (S 2 ) x {1} and projected centrally on the unit upper Hemisphere H. Its equator £ is the 
sphere at infinity. 



hemisphere H is not compact. Following [19], the Poincare compactification induces a com- 
pactified equation on % as well as an equation on the sphere £ at infinity, which is the one 
we are interested in. 

We begin from the equation (11.4)) . which can be simplified in our simple settings by using 
u i — y r~ x l 2 u to 

(7.38) d r u = u 2 Au + ^u 3 . 

Sec [32J for details. Then the Poincare compactification induces after normalization the 
following equation on the sphere at infinity: 

(7.39) a,X^Ax+^-(x°Ax + ^,x)x, IMI=1, X ^, 

where (., .) is the scalar product on L 2 and the time variable r corresponds to a normalization 
preventing the trajectories from hitting the sphere at infinity in finite time. The previous 
equation (17.391) is not a full PDE because it contains nonlocal terms coming from the scalar 
product. Now we can interpret the self-similarly rescaled equation (II . 14[) ( which differs from 
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fll.ISp only by a multiplication by a scalar and a translation of the origin) as an equation on 
the sphere at infinity. Let v be a L 2 (S' 2 )-bounded solution of equation (j!.14p 



Then \ = i^fL G £ solves the equation on the sphere at infinity (17.391) up to a factor that 



The heteroclinic connections discussed previously are objects in the sphere £ at infinity, and 
influence the blow up behavior of our original Problem as depicted in Figure I7.3[ which is a 
3-d caricature of this infinite dimensional hemisphere (ball) %. Here we see that trajectories 
with finite initial conditions in the interior of the Poincare hemisphere T-L are caught by 
the heteroclinic connection between the trivial equilibrium at infinity and a manifold Tv* of 
anisotropic equilibria at infinity, and hence explode anisotropically. 

Now that we have sketched the dynamics at infinity of the original Equation (jl.4p . let us 
interpret its geometrical meaning for the manifold M at its upper boundary {r{\ x S 2 . In 
fact, the blow-up behavior reflects a degeneracy of the area-radius r, but not a degeneracy 
of the metric g of the manifold M at its upper boundary {ri} x S 2 . Let us recall that the 
metric on the manifold M is given by 



where 00 is the standard metric of the round sphere. Now we introduce a coordinates s 
proportional to the geodesic length as 



d t v = v 2 l\v + -v 3 - v. 




(7.41) 



2j 2 1 2 

g = u dr + r oj, 



(7.42) 




Hence 



(7.43) 



dr = (1/r 



l)ads, 



and the metric belonging to the (asymptotically) self-similar blow-up u 



(f (Vr-1 



)) *„ 



reads in these coordinates: 
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Figure 7.3: This Picture caricatures the infinite dimensional Poincare hemisphere H of 
L 2 (S 2 ) in 3-d. The single point denoted by 1 is the isotropic equilibrium at infinity. The 
ellipse denoted by Yv* stands for a manifold of anisotropic equilibria. The straight dashed 
lines coming from the interior of the Poincare hemisphere are self-similar blow up solutions 
with profile 1 and v* respectively. The sphere £ at infinity contains a heteroclinic 1 — >■ v* 
that influences the blow up behavior of non self-similar blow up solutions. 

In this form, the metric is continuous and nondegenerate up to and including the upper 
boundary of the manifold M: g £ C°°(M) D C°(M), no matter which profile for v is used. 
It may be the trivial equilibrium of the self-similarly rescaled Equation (11.181) . one of the 
equilibria or heteroclinics we exhibited nearby via equivariant bifurcation theory, or any 
solution converging to any other equilibrium v of the self-similarly rescaled Equation ( 11. 18ft . 
The forward time evolution of a heteroclinic connection between two equilibria v\, v 2 of 
Equation (I1.18P has just been discussed. Let us now discuss the backward time evolution 
t — > — oo alias r — > — oo through r = 1 — exp(— 2t) of a heteroclinic trajectory v(t) from 
equilibrium v\ to equilibrium v 2 . Replacing the latter change of time/radius variable by f = 
1 — exp(— 29) exp(— 2t) corresponds to a time shift t t— > t+6. Let tq be the initial radius of our 
construction. The first change of radius/time variable puts ro at time to = l n ((1 — r o) _1 ^ 2 ) ; 
while the second puts ro at to + in the second case, v(to + 9) determines the shape of the 
metric at initial radius r : v{t Q + 9) comes arbitrarily near equilibrium v\ as 9 approaches 
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— oo. Once 9 is fixed, the equilibrium v 2 at the other end of the heteroclinic gives the shape 
of the metric at the blow up radius r% = 1 t,t + 9 — > oo for any fixed 9. 
Let us now discuss the question of the center of M, that is, the limit where the lower 
boundary disappears as r — > 0. Because we chose to study bifurcation with respect to the 
parameter A, this imposes the prescribed scalar curvature to be R(r) = ^4^. Therefore, the 
scalar curvature R, and with it at least one sectional curvature, concentrates at r = 0: we 
cannot expect M to be regular at the center. 

This is seen again by considering the coordinate s proportional to the geodesic length, and 
for simplicity the isotropic self-similar blow-up profile v = 1, with blow up radius r% = 1. 
Near r = 0, the radial component of the metric can be approximated by j^- ~ r, so that 
the metric and the geodesic length are given by 

g = rdr 2 + br 2 u 

(7.45) s = y/fdf = |r5 

g = ds 2 + bssoj 

The singularity manifests itself in the term si: the discrepancy between the exponent | and 
2 excludes regularity of the center. 

Our motivation was to construct M as a piece of initial conditions for the Einstein equations. 
In view of this, curvature singularity at the center may not be desirable. To avoid such a 
situation, one may want to fix a r > and use a different method for prescribing the scalar 
curvature for r < tq. 

Even if the isotropic scalar curvature R(r) = which we have prescribed in this work for 
simplicity in the arguments involving bifurcation theory is rather primitive, we are still able 
to show that a rich anisotropic and only asymptotically self-similar behavior develops. The 
isotropy of the metric component in radial direction can break in several ways as the bound- 
ary {ri} x S 2 of the manifold M is approached, depending on the choice of the bifurcation 
parameter A. This method may be applied to more realistic situations where u and r 2 R do 
depend on r. In this case, it may be possible to apply the equivariant bifurcation theory 
arguments on a thin collar [ro, ri) x 5 2 region, and complete the manifold by prescribing a 
scalar curvature leading to a regular center in the center region [0, ro] x S 2 . 

Let us briefly address the behavior of the metric g at r = 0, according to our construction. 
The parabolic scalar curvature equation 11.41 is equivariant, separately, with respect to each 
of the involutions rf> -r 6 [-1,1] and ip -H- — y? G X = S 2 . Our transformation u i— > v to 
the self-similar rescaled equation (11.131) preserves both involutions if we define 

(7.46) v = y/l/\r\ - 1 u, 
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to accomodate negative r. A simple reflection 



(7.47) z/(— r,p) := v(r,p), 

for r > then provides a consistent extension to negative radial arguments. In terms of 
smooth elements v in a strong stable manifold W ss of fll . 18j) we get 

(7.48) u = (l/\r\ - \y ll2 v = |r| 1/2 (l + -|r| + . . >. 

With z/ = (A/2) [v + 1), and u regular in t — r + . . . at r = we obtain an expansion 

(7.49) u 2 = \r\(l + \r\ + ...)?- (1+v) 2 . 

A 

This identifies the metric 

(7.50) g = u 2 dr 2 + r 2 u 

to possess a quadratic cusp at r = 0. Nonetheless, Definition (17.471) provides a valid extension 
of u beyond r = into negative r. 

Suppose we then insist, in addition, that (r, p) and (— r, — p) denote the same point on the 
3- manifold M = [0,ri] x S 2 as common in polar coordinates at the point r = where the 
cusp singularity resides. Then we have to require antipodal symmetry 

(7.51) z/(r, — p) = v(—r,p) = is(r,p) 

for v and consequently for v. This requirement amounts to isotropy of v under —id G 0(3), 
in the language of Section [2j 

(7.52) r„ > 7L\ = (-id). 

By our 0(3) bifurcation analysis of Section 3, this requirement is automatically satisfied for 
all symmetry breaking branches (A(s),t>(s)) which emanate at A(0) = A^ = l[l + 1) with 
even I = 0, 2, 4, . . .. But (I7.5ip .( 17.52p is just as automatically violated at odd t. From this 
point of view, even I are preferable because they do not need any further treatments for 
< r < r\ and because the simple reflection (17.471) defines a consistent, and symmetric 
extension to — r\ < r < 0. 

As noted previously, the initial data for the Einstein equations have to be completed by a) a 
regular center, and b) a region matching the blow up boundary {rx} x S 2 with the apparent 
horizon of the black hole. The latter can be constructed with anisotropy, for more details see 
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Both sides of the matching region, {ri} x S 2 and the apparent horizon are critical 
points of the area functional, the latter being a minimum of Morse index 0, and the former 
a critical point of index greater than 1. To see this, let us notice that at blow up radius 
r = 1 , the mean curvature H = H/ru oc (1 — r) 1 / 2 /r 3//2 vanishes, and this implies that the 
sphere S 2 x 1 is a critical point of the area functional A, i.e. a so called minimal surface. If 
we denote by F T : E = S 2 — > M a smooth family of spheres parametrized by a parameter r, 
such that d T F is a multiple of the outer normal N to S T = F T (S 2 ) with factor ip, d T F = ipN, 
then (see HO]) 



dA f 

(7.53) = - / ipHdA. 

dr J S 2 

To analyze the type of critical point we are facing, we need the second variation of the area 
functional at r = 0, which reads in our simple settings 

d 2 A 



(7.54) 



/ |V^| 2 + V 2 (l-R/2). 
Js 2 



The self adjoint operator associated to this functional is A — (1 — R/2) = A + A/2, and the 
stability index of a critical point of A is determined by the number of positive eigenvalues 
(counted with their multiplicity) of this operator. Hence we get for A G (2A^, 2A^ + i) that the 
stability index of the minimal surface at blow up radius is = 2£(£ + 1) + 1 
Once the completions a) and b) have been performed, initial conditions with anisotropic 
metrics will be provided for the Einstein equations. We will then be in a position to ask 
about their time evolution and study, in particular, how the anisotropy affects the long time 
behavior. 
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